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I. INTRODUCTION 

Recent numerical studies jl|] have provided strong evidence that /7(l)-symmetric, vacuum spacetimes on x R 
generically develop Mixmaster- like, oscillatory singularities of the type predicted long ago by Belinskii, Khalatnikov 
and Lifschitz (BKL) These results confirm numerically some of the most surprising features of the BKL 

prediction, namely that nearby spatial points are effectively decoupled in their asymptotic metric evolution and that 
the metric variables at each of these points evolve, at least qualitatively, like those of a Mixmaster spacetime. 

Several years ago B. Grubisic and one of us (V.M.) made an analytical effort to generate some exact vacuum 
spacetimes which were spatially inhomogeneous and which were expected to exhibit the sort of oscillatory singularities 
which have since been seen in the numerical studies [|l|. That effort was not completed at the time since it was not 
realized that several seemingly intractable integrals actually cancel in the course of the calculations leaving only 
elementary computations to be done. We shall therefore complete that project here and use the results to compare, 
in a more quantitative way, the numerical results with some exact oscillatory singularities. 

To generate new solutions having Mixmaster-like oscillations, we begin with the actual Mixmaster solutions and 
apply a standard solution generating technique. We choose one of the Killing fields shared by the Mixmaster family 
and treat it as the generator of a spacelike U{1) action on x R, ignoring the presence of the other Killing symmetries. 
We compute the twist potential associated with the chosen Killing field and reexpress the field equations, in a well- 
known way as a Kaluza-Klein reduced system on the base manifold 5*^ x R of the 5'^-bundle 6*^ x R ^ 5^ x R. 
The field equations on the base take the form of 2 + 1 Einstein gravity coupled to a wave map whose target space is 
the hyperbolic plane, H^. The isometry group of this latter space, 5'L(2,R), acts on the base fields in a natural way 
so as to transform the given solution to a family of potentially inequivalent solutions. 

By a careful choice of the applied group element one can arrange that the transformed solution either lifts to the same 
bundle defined for the original spacetime or perhaps to a different one (e.g., the trivial bundle, 5^ x 5^ x R ^ 5^ x R, 
or a "squashed sphere", x R ^ 5^ x R). Typically, the new solutions will preserve only the Killing field that 

generates the common U{1) action and not preserve those Killing fields of the seed solutions which fail to commute 
with the chosen U{1) generator. Thus the new solutions are expected to be spatially inhomogeneous and yet to 
exhibit Mixmaster-like oscillations inasmuch as their metrics are parametrized by the same functions appearing in the 
Mixmaster seed metrics themselves. 

A previous application of this technique involved transforming an infinite dimensional family of "generalized Taub- 
NUT" spacetimes defined on x R, which have smooth Cauchy horizons at their "singular" boundaries, to a new 
family of curvature singular spacetimes defined on 5*^ x 5*^ x R ||] . Because of the special nature of the seed solutions in 
this case, the transformed solutions developed only velocity dominated singularities and never exhibited Mixmaster- 
like oscillations. A new technique based upon expressing the Einstein evolution equations in a so-called Fuchsian 
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form seems capable of significantly enlarging this set of rigorous, [/(l)-symmetric, curvature singular cosmological 
spacetimes but, so far, is also only capable of yielding velocity dominated singularities |l^,|ll|]. So far as we know the 
solutions presented for the first time here are the only known exact inhomogeneous vacuum spacetimes which exhibit 
Mixmaster oscillations. Though only a finite dimensional family they presumably display behavior representative of 
more general, /7(l)-symmetric vacuum spacetimes and thus warrant comparison with numerically produced U{1)- 
solutions. Making such a comparison is the second main aim, after producing the solutions themselves, of this paper. 
As a byproduct of this work, we also resolve a potential paradox that was pointed out in Ref. 0. There it was 
shown that every /7(l)-symmetric vacuum spacetime admits a certain gauge invariant conserved quantity which is 
expressible purely locally in terms of the instantaneous Cauchy data for that solution and serves as a Casimir invariant 
for the 5'I/(2,R) action. For generic U{1) solutions this quantity is known to be non- trivial but, if non-trivial for the 
Mixmaster subfamily, would seem to contradict the anticipated "chaos" of the Mixmaster dynamics [|l2|-|l^ . The only 
sensible resolution, as was discussed in Ref. 0, is that the quantity actually vanishes on the Mixmaster subfamily. 
This we find to be the case by explicit calculation. 

The inhomogeneity in our transformed solutions is produced, roughly speaking, by the fact that we choose to reduce 
with respect to a Killing field which fails to commute with the remaining Killing fields of the seed metric. This is 
unavoidable with the generic Mixmaster solution but special cases such as the Taub-NUT metrics allow for different 
possibilities. The additional Killing field admitted by Taub space commutes with all the generators and is preserved 
upon reduction with respect to one of these (non-abelian) generators. The resulting spacetime has therefore (at least) 
two commuting Killing fields and is thus a special case of the so-called Gowdy family of spacetimes. By contrast 
one could instead choose to reduce with respect to the additional, commutative Killling field but, in this case, all 
the symmetries are preserved and one arrives, as was first shown by Geroch at only the Kantowski-Sachs (i.,e. 
locally interior Schwarzschild) spacetime. 

One might wonder if the "new" solutions we produce are really inhomogeneous at all or perhaps because of their 
expression in an unusual gauge, are merely homogeneous solutions in disguise. We shall use the Gowdy transform 
of Taub space mentioned above, to show that this is not the case — the new solutions are not in general globally 
homogeneous. 



II. MIXMASTER SPACETIMES 



The Mixmaster spacetimes are spatially homogeneous vacuum metrics on /S^ x R whose line elements can be written 

ds^ = -N\t)dt^ + A\t){&y + B\t){a^)^ + C\t){a^)^ (1) 

Here the {a*} are a global, analytic basis of one- forms on expressible in terms of the usual Euler angle coordinates 
{x\x\x^} = {e, v>, V'} e {[0, tt), [0, 2n), [0, 47r)} by 

= cos (f dO -\- sin sin (p dip^ 
(j^ = — sin (f dO -\- sin cos (p dip^ 

= dif^ cos dip . (2) 

These forms, and therefore the above line element, are invariant with respect to the U{1) action on generated by 
the Killing field Xs = ^ whose orbits yield a Hopf fibration of 5^, i.e. make into a principal fiber bundle over S"^ 
with bundle projection given by 

7r^:S^^SMe,ip,xP)^{e,^). (3) 
Of course the Mixmaster metrics are invariant with respect to a full SU (2) action generated by Killing fields 

d d d 

Xi — cosip—- + CSC Osimp— cot Osimp—-^ 

oO Oif dip 

X2 = — sin?/^— - + CSC cosip— cot 6 cosip——^ 

06 0(f dip 

^3 = ^ (4) 

but, for the transformations we shall consider, only invariance with respect to X3 = ^ will in general be preserved. 
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The equations of motion for the Mixmaster solutions are most simply expressed in a gauge for which N = ABC 
where they take the form 



(5) 

and are to be supplemented by the Hamiltonian constraint 

AB AC BC 

^B^^C^- 2{A^B^ + B^C^ + A^C^)] ^ 0. (6) 

In terms of the Misner anistropy variables 

Q ^ ga+/3+-V3 (3- ^ 

C = e"-2^+, (7) 

and the chosen gauge condition is A/" = e^". We now rewrite the line element in the /7(l)-symmetric form developed 
in Refs. ||^ and |l^]. Taking {x^} = (f} and noting that the shift vector vanishes we express ds^ in the form 

ds^ = e-'^^{-N'^df + gabdx^'dx^} 

-^e^^id^lj + cos i9 di^ + Mx^'Y (8) 
where e^^ is the scalar field ^ • ^ given explicitly by 

e^^ = A^ sin^ 6 sin^ + 5^ sin^ l9 cos^ if 

cos^ 6> . (9) 

Since 7 is invariant with respect to the U{1) action generated by X3 = ^ it induces a function on the quotient 
manifold x R//7(l) S*^ x R which (with a slight abuse of notation) we shall also designate by 7. In a similar 
way one finds induced upon the quotient manifold 5^ x R a Lorentzian metric 

da'^ = -N'^dt^ + Qahdx^'dx^ (10) 

and a one-form field 

(3 = Mx"" (11) 

where {a, 6, . . .} = {1,2}. These forms are slightly specialized because of the vanishing of the shift vector field in ds^ . 
The most general /7(l)-symmetric line element would yield 

dcr^ = -N'^dt^ + ~gah{dx'' + N''dt){dx^ + N^dt) 
p = padx"" + podt. (12) 

The explicit formulas for the 2 + 1 Lorentzian metric da^ and the one-form potential /3 may be read off upon 
expressing ds'^ in the form of Eq. (|^). One finds that 

TV = TVe^, 

gee = cos^ 6{A'^ cos^ (p ^ B'^ sin^ cp) 

^A^B'^ sin^ i9, 
g^^ = sin^ 6>(A^ sin^ (p ^ B'^ cos^ cp), 
geif = —C^ {A^ — B'^) cos (p sin (p cos 6> sin 6>, 
(yA^ — B^ ) COS p) sin sin ^ 



e 
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and computes, for example, that 



Vdet (2)^ = ABCsinOe^. (14) 



As in Refs. [|| and we introduce the momenta {p^e^^n^^} conjugate to {j^ (3a^ Qab}^ which, taken together, 
parameterize the fuh 3 + 1 spatial metric gij and its conjugate momentum tt*-^ . For the case of vanishing shift the 
formulas relating the momentum variables {p, e^,7r^^} to the metric variables {j^Pa^9ab} are given by 
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For the Mixmaster metrics one computes that 

2 sin^ 

= — — — {BCA t sin (p cos (p — ACB t sin ip cos (/:?}, 

= '^^^^^lpl{ABC^t - BCA^t sin' ^ - ACB^t cos' (^}, 

p = ^^f,^,^'""^ [AA t sin' sin' ^ ^ BB t sin' ^ cos' ^^CCt cos' ^] . (16) 

These momenta (along with tt^^ which we shall not need explicitly) project to yield smooth tensor densities on the 
base manifold 5'^ x R and one easily verifies that e\ = which is one of the components of the 3 + 1 momentum 
constraint. 

Note that the U{1) connection one- form on x R given by 

A := cD^ + /? = # + cos Odif + fSadx"" (17) 

does not project to yield a one-form on the base but that the difference between this connection one form and the 
reference one form uj^ := dip -\- cos Od(p does yield a one-form (namely p = padx^) which projects to the base. Even 

though A itself does not project to the base, its exterior derivative d\ (i.e., the curvature of the connection A) does 

project. Pulling back the induced two- form to a t = constant slice of the base manifold and computing its dual, one 
gets a scalar density f defined by 

= G"' I3afi + sin e (18) 

whose explicit form is 

^_{A^- 5') ^.^3 ^[^2 ^^g2 ^ _ ^2 g-^2 



(e27)2 

sin (9 



[A' + 5' - C' + sin' e{C^ - A' cos' - 5' sin' (^)]. (19) 



(e27)2 

One computes on an arbitrary t = constant slice of the base manifold, that 



B := f = 47r. (20) 

The value 47r reflects the particular bundle 5^ x R ^ 5'^ x R under study and would be the same for any U{1)- 
symmeteric metric defined on this bundle. 
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Taking into account the equation e^^ = satisfied by and the fact that S*^ x R admits no non-trivial harmonic 
one forms, we now introduce the "twist potential" function co (a salar field on 5*^ x R) by imposing 



N 

These equations are self consistent and yield the solution 

.2 /) 



^^e-^,. (21) 



-{-ABC^t + BCA^t sin^ ^ + ACB^t cos^ ^} 



N 

-^k{t) (22) 
which is unique up to the additive constant ko := k{to) where k{t) is the function defined by 

m = k{to) + l^'dt' (^-^) C\A' + C'). (23) 

As discussed in Refs. [|| and the fields {j^u^da'^ = —N'^dt^ + gabdx^dx^} induced upon the base manifold 
5*^ X R satisfy a 2 + 1 dimensional system of Einstein-wave map equations for which the target space of the wave 
map is hyperbolic two-space (endowed with global coordinates {j^oo} and the natural metric dh^ = 4dj^ + e~^^duj'^). 
As a consequence of the 5'L(2,R) isometry group of this target space the Einstein-wave map system admits three 
independent constants of the motion which serve as the Hamiltonian generators of the action of SL{2^ R) on the phase 
space of fields {7,^, cj, f, ^^65 tt^^}. These conserved quantities are given explicitly by the integrals 

i:= / (2cjf + p), 



B:= f, 

C := [ (rV^-c.2)-pc.), (24) 

and we have already noticed that B = Att for the Mixmaster spacetimes in particular. In fact B would take this same 
value for any /7(l)-symmetric vacuum metric on S*^ x R but for other bundles over the same base the value would 



(as discussed in Refs. [|[ and |16|) be modified to B = Ann where n is an integer determining the Chern class of the 
bundle. In particular, for solutions on the trivial bundle 5*^ x S*^ x R, n would vanish whereas if n = 2, 3, . . . , the 
bundle would correspond to various "squashed spheres" rather that a true . 

Note that, in view of the integral expression for k(t) arising in the formula for cj, both A and C are non-local 
in time. The same feature occurs in more general U(l) symmetric solutions but this non-locality cancels from the 
Casimir invariant 

k-=A^^4.BC (25) 

which, however, vanishes identically for the Mixmaster family of solutions (though not in general). The vanishing 
of K resolves a potential mystery pointed out in Ref. whereby a non-vanishing, local, constant of the motion for 
Mixmaster metrics would seem to contradict their empirically observed "chaotic" properties. 
More specifically one finds, for the Mixmaster metrics, that 



fABC\ [At Bt 



B = 47r, 

so that K = 0. The non- locality of A and C sidesteps any conflict with the observed "chaos" in Mixmaster solutions 
since, in fact, any Hamiltonian system will admit such non-local constants of the motion. To see this (even for a 
chaotic system) simply time integrate Hamilton's equations and express the initial values of the canonical variables 
in terms of time integrals of their driving "forces." 
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III. THE NEW SOLUTIONS 



To generate new solutions of Einstein's equations from a given one (such as a Mixmaster solution) we choose an 
element g G 6'L(2,R), 

9=(yl J), ad-hc=l (27) 
and transform the fields {7,cc;,p, f} according to 

p27, _ ^ 



[C2(cj2+e47)+2cdcj + d2]' 

_ ac(cj^ + e^^) + {ad + hc)uj + hd 
. _ {p[c^(cj^ - e^^) + 2cdcj + d^] - r[4e^^{cd + cc;c^)]} 

= pic^co + cd) + r[d'^ + c^(cj2 - e"^^) + 2cdu;], (28) 

while leaving {^^6, tt'^^, TV, 7V^} invariant. The induced transformation of the conserved quantities A^B^C (by the 
so-called co-adjoint action of 6'L(2,R))is found to be 

Ag = {ad + be) A + 2bdB - 2acC, 
Bg = d^B - c^C + cdA, 
Cg = a^C - \?B - abA, 

kg=k={Agf^4BgCg. (29) 

To avoid a trivial transformation we shall require that c be non-zero and, to ensure that the transformed solution 
lifts to an bundle overs'^ x R, we shall demand that 



Defining 



we see from Eq. (|27| ) that 



Bg =47rn,n = 0,1,2,... . (30) 



A = S7rk{t) + S7ri{t) = S7rk{to) + S7ri{to) = 87r(A:o + 4), 
B = 47r, 

C= -4^(^0+4)'./ (32) 
Setting Bg = 47rn, n = 0, 1, 2, . . . gives the restriction 

[d^c{ko^io)]^ =n>0 (33) 

or, equivalently, 

d + c(A:o+4) = ±n^/^ (34) 

which can always be solved for ko since c 7^ 0. 

Exploiting the fact that A, hence also k{t) + i{t)^ is conserved one finds that the integral occurring in the formula 
for k{t) can be expressed as 



ABc) 



C^{A^ + B^-C^) = -{e{t)-eo) (35) 
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(this is also easily verified upon differentiation by using the equations of motion (||)). Using this result, one can easily 
show that 

+c [-ABC^t + BCA^t sin^ ^ + ACB^t cos^ ip] 



ABC ( A^t Bt 



N \ A B 



(36) 



With this and Eq. (^ for e^^ one easily evaluates the transformed field variables {e'^^s ^cjg^Pg, Tg} using Eq. (p8|). The 
new spacetime metric thus takes the form 

dsl = e-'^^^{-N'^dt^ + gabdx^'dx^} 

-^e^^^idijj -^n cos Odif + (3i^g)adx''}'^ (37) 

where however, P(g)a remains to be computed. As discussed in Ref. [||, P(g)a can be expanded (via the Hodge 
decomposition for a one- form on 5^) as 



(38) 



where fig = and r] and 5 are suitable functions defined on S'^ . The equation for (^{g)a dx^ is 

fg =G"^/3(^),,5 + nsin ^ (39) 

which, upon substitution of the decomposition (|38|), becomes 

f^-n sin = {flgg''^r]^a),b (40) 

a Poisson equation for 77 for which the necessary and sufficient integrability condition is ensured by Eq. (|3Q|). This 
uniquely determines 7^, at fixed up to an arbitrary additive constant and leaves 6 arbitrary. The presence of 6 
reflects the freedom to make an arbitrary coordinate transformation of the form iIj ^ ip -\- S without affecting the 
[/(l)-form of the spacetime metric. 

The time development of P(g)a can now be obtained by integrating the (zero shift) evolution equation 



9) 



e-^^^ -gate'' ug^c (41) 



with Jg^ujg determined as above. 

Equations (|39| ) and ( p]) are consistent with each other by virtue of the Hamilton equations satisfied by Vg^uOg. Note 
that whereas we have used the actual metric gab in defining a Hodge decomposition of I3(^g)a-, any smooth metric on 
S'^ could have been used instead. Furthermore one could have used Eq. (ji^) to determine r] at an arbitrary time and 
then adjusted the time dependence of 3 to impose the zero shift condition which is implicit in Eq. (|4l|). In either case 



the new metric (37) will satisfy the vacuum field equations on the chosen bundle over S'^ x R. 

One might still wonder how we know that the transformed solutions are genuinely inhomogeneous. Could they not be 
merely homogeneous solutions disguised through the choice of a time slicing that is not adapted to the (hypothetical) 
homogeneity? To show that this is not the case, in general,we shall examine a special case for which the transformed 
solution has a hypersurface of time symmetry at t = to, i.e., has k[^^ \t=to = 0. To arrange this, we choose the seed 
solution to have this property and make a careful choice of transformation parameters so that the desired feature is 
not destroyed by the 5^(2, R) transformation. We then show that the transformed spatial metric g\^^ \t=to is not 
homogeneous as it would have to be for the resulting spacetime to have this property. The key point here is the fact 
that on any compact slice, having constant mean curvature the first and second fundamental forms {g\^\K\^^} \t=to 
would both have to be homogeneous in order that the spacetime have this property. 
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Consider a Mixmaster solution for which A{to) = B{to) = C{to) = 0. This spacetime has the t = to shce as a 
surface of time symmetry and, because of the Hamiltonian constraint, must satisfy 



C{to) = ±{A{to)±B{to)), 

A(to),5(to),C(to) >0. (42) 

To maintain this property we choose the trivial target bundle S'^ x x H hy taking n = 0. We further simplify 
the computations by choosing A{to) = B{to) > and C{to) = 2A{to) and find that the transformed metric at t = to 
satisfies 



gee 


It= 




geip 


l*= 


=to 


gipcp 


\t- 






\t-. 


=to 


e27« 


\t- 





- J4/ 



1 COS 



= ^0 = 

1 



cM2(to)(sin2 6» + 4cos2 6') " 
Thus the new spatial metric induced a.t t = to on S"^ x 



■4A'^(to)cos6'sin2 0dip, 

(43) 



difg-^ \t=to = c2A2(to)(sin2 + 4cos2 0){4A'^ito) sin^ ddip^ 
+A'^{to){sm^ 6 + 4 cos^ e)de'^} 



-— -— — i— — — -{d^ - ic'AHh) cos^sin^ 0d^}' 
c^A^(to)(sm + 4cos^ 0) 



(44) 

A straightforward computation of RijR^^ (the square of the Ricci tensor of this metric) proves that the resulting 
spacetime is not homogeneous. Indeed, the only vacuum homogeneous solution on S'^ x x H is known to be the 
K ant owski- Sachs universe which does not have a hypersurface of time symmetry. It is possible to get the Kantowski- 
Sachs metric upon transformation of a Taub metric (i.e., a Mixmaster solution having A{t) = B(t)) but to do so one 
must reduce with respect to the "extra" Killing field the Taub metric possesses, rather than with respect to the 

common Killing field ^ of the Mixmaster family as we have done. The extra Killing field ^ commutes with all the 
Killing symmetries of the Taub solution and allows all of these symmetries to be preserved upon reduction [ p!5| . 

The Taub metric used in the example above, is known explicitly but exhibits no BKL type oscillations. To see 
such oscillations in an inhomogeneous setting, we combine a previously developed code for solving the Mixmaster 
equations of motion with the transformations discussed above. Our results are discussed in the following section and 
compared with results derived from a general [/(l)-syummetric, vacuum Einstein code. 



IV. NUMERICAL RESULTS 



Elsewhere we have shown that even the homogeneous Mixmaster model reproduces the local behavior seen in generic 
/7(l)-symmetric cosmologies Jl^. From Eq. (^), it is clear that 7 is dominated by the largest of the Mixmaster scale 
factors ^4, 5, or C. The local oscillations seen in 7 in the /7(l)-symmetric models are interpreted as follows: Assume 
that the BKL approximate description of a homogeneous Mixmaster model as a sequence of Kasner epochs is valid. 
In a given approximate Kasner epoch assume A > B > C and that A is increasing. Then 7,^ > for A^t > while 
B^t and C,t are less than zero. The usual Mixmaster bounce changes the sign of A^t and thus of 7,t. However, after 
the bounce, either B^t (within an era) or C,t (at the end of an era) becomes positive. When the growing scale factor 
surpasses the decreasing A, 7,^ will start to grow again since it will now track the new dominant scale factor. A 
similar analysis indicates that the remaining "dynamical" variables, cj, and f, depend on an order unity ratio of 
scale factors and thus do not oscillate, as 7 does, between order unity and exponentially small values. (Here we shall 
use "order unity" to mean some finite value which is not exponentially small.) 

In our previous numerical simulations of generic /7(l)-symmetric cosmologies on x R, we noted that the os- 
cillations in 7 could be interpreted as bounces off the potentials Vi = ^f^e^^ and V2 = g^^e~^^u,a^ib' For a 
Mixmaster solution, Vi is exponentially small unless e^^ is of order unity while V2 is exponentially small unless the 
two largest scale factors are approximately equal to each other [|l^. This is clearly consistent with a presumption 
that the generic models exhibit local Mixmaster dynamics. 
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To explore the nature of the new inhomogeneous f/(l)-symmetric models, we note that the transformed variables 
uOg^ Pg^ and will remain of order unity (i.e. they will not oscillate between exponentially small and order unity 
values) because the right hand sides of Eqs. ([2^b)-(p8|d) are always of order unity. On the other hand, 7^ is dominated 
by the behavior of the oscillatory 7 since the denominator on the right hand side of Eq. ([2^a) is always order unity 
while the numerator oscillates. 

To explore the differences between our new solution and the Mixmaster seed solution, we construct the new solutions 
as follows: First use the algorithm of Berger et al to obtain a numericallly generated Mixmaster model. This 
code is known to solve the Mixmaster ODE's with machine-level precision and can follow hundreds of bounces. The 
presumed stochastic properties of such a model imply that almost any Mixmaster initial conditions will yield generic 
Mixmaster behavior. Thus, we need only consider a single Mixmaster trajectory. Next, Eqs. (||), ([l6|), (|l9|), (|2^), 
(p3|), ([31]), and ( |35| ) are used to numerically evaluate 7, cj, p, and r from the numerically generated sequence of values 
of the BKL scale factors and their time derivatives. Finally, for a representative choice of the /SL(2,R) parameters 
and, e.g., n = 1, the transformed variables 7^, o;^, Pg, and Tg are computed using Eqs. (p8[). 

In Figures 1-3, we compare the Mixmaster and transformed 7 and cc; at a representative spatial point for typical 
Mixmaster seed and set of /SL(2,R) parameters. Note that, in Fig. 1, the original and transformed 7's become 
indistinguishable after only a small number of Mixmaster epochs. It is clear that this will be so from Eq. (^) for 7 
and Eq. (^sja) for the transformation. Since 7 and 7^ are found from the logarithm of Eqs. and (|2qa), both 7 
and jg will be approximately equal to the logarithm of the largest scale factor and depend only logarithmically on 
the spatially dependent function associated with it. On a finer scale, in Fig. 2, the difference between the solutions 
(especially near the "bounce" where j ^ jg ^ 0) may be seen. On the other hand, as is seen in Fig. 3, uj and Ug are 
always of order unity and may easily be distinguished. 

Figure 4 demonstrates the close link between Mixmaster dynamics and the oscillatory behavior observed in our 
studies of generic [/(l)-symmetric models and should be compared to Figs. 2-6 in [||. It shows the oscillations of 
7p (or essentially equivalently 7), Vi, and V2 at a representative spatial point — reproducing the behavior seen in our 
simulations of generic /7(l)-symmetric cosmologies. Since we know that these oscillations indicate local Mixmaster 
behavior in the new solutions, we can infer that the observed oscillations in the generic models also indicate local 
Mixmaster dynamics. 

Since 7 is the key variable in the /7(l)-symmetric models and 7 ^ 7^, one may then ask where these new U{1)- 
symmetric models differ from both Mixmaster and generic /7(l)-symmetric models. First, we emphasize that, except 
at special values of the spatial coordinate angles, there are no qualitative differences attributable to spatial topology. 
The Mixmaster spatial dependence of course represents a realization of the Bianchi type IX symmetry. From Eq. (||), 
it is clear that three distinct spatial patterns will appear in 7 (in the logarithm) depending on which scale factor 
dominates. In Fig. 5, we compare the spatial dependence of 7 and 7^ for 12 epochs of the seed Mixmaster solution. 
The epochs are arranged according to the dominant scale factor. The numerical scale in each frame is chosen so that 
the average value of 7 or 7^ is the centroid. (If this were not done, no spatial dependence would be visible.) From 
Eqs. (^ and (22), 7 and co have three possible spatial dependences. The SL(2^ R) transformation of Eqs. ( p8| ) clearly 
mixes the spatial dependence of 7 and uo to form 7^ and uOg. In Fig. 5, we see the evolving spatial dependence of 7^. 
This is additional evidence that the new solutions are spatially inhomogeneous. 

In generic /7(l)-symmetric models, one could qualitatively interpret the asymptotic approach to the singularity as 
the evolution of a different Mixmaster model at every spatial point. In particular, the Mixmaster epochs have spatially 
dependent durations — bounces at different spatial points occur at different times. In contrast, our new solution is 
characterized as is Mixmaster itself by spatially independent epoch durations since the new solution bounces only 
when the seed solution does so. While one could modify the spacetime slicing to yield spatially dependent epoch 
durations, one would expect to be able to detect the difference between between a single underlying Mixmaster seed 
in the new solutions and a continuum of approximate Mixmaster solutions in the generic case. 
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FIGURE CAPTIONS 



Figure 1. Comparison of 7 and 7^ at a typical spatial point. The Mixmaster seed solution has initial values 
/3+ = -0.9847899998176387, p- = 0.09987655443789, Q = -8.00000000000000, /3+ = -3.632980009876544, = 
4.58987654433567878, and the Hamiltonian constraint solved for (l. The 51/(2, R) parameters are a = 1, 6 = 1, 
c= 10000, and(i= 10001. 

Figure 2. Detail of the comparison of 7 and 7^. To emphasize the approach of 7^ to 7, data from later in the 
simulation of Fig. 1 are shown. The actual, saved data values are indicated by the x and + symbols. 

Figure 3. Comparison of co and cOg for the same models as in Fig. 1. Note that cOg appears to decrease to zero. This 
is due to the fact that choice of 5^(2, R) parameters causes \ug\ ^ 10~^ if |cc;| ^ 1. 

Figure 4. New solution as an inhomogeneous /7(l)-symmetric cosmology. As in previous studies of generic U{1)- 
symmetric cosmologies, 7^, Vi, and V2 are shown at a typical spatial point. 

Figure 5. Evidence for the spatial inhomogeneity of the new solutions. The spatial dependence of 7 and 7^ is shown 
for the (cos^, cp) plane in a series of side-by-side frames arranged in three separate panels. Each pair of frames shows 
the spatial dependence of 7 and 7^ respectively during an approximate Kasner epoch of the seed Mixmaster solution. 
The panels are grouped according to the identity of the dominant scale factor in the spatially homogeneous solution 
rather than sequentially. According to Eq. (||), 7 will have the spatial dependence ln(sin6> sin(/?), In(sin6> cos (/:?), or 
In(cos^) depending on whether A, 5, or C respectively is dominant. In each of the three panels, the 4 left-hand 
frames reproduce one of these three spatial dependences with, reading from left to right, 5, C, or A dominant. In 
each case, the accompanying right-hand frame represents the spatial dependence of the corresponding 7^ for that 
epoch. In every case, the numerical scales for 7 and 7^ have been centered on their average values to enhance the 
visibility of the spatial dependence. 
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